Model inspired by population genetics to study fragmentation of brittle plates 
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We use a model whose rules were inspired by population genetics, the random capability growth 
model, to describe the statistical details observed in experiments of fragmentation of brittle platelike 
objects, and in particular the existence of (i) composite scaling laws, (ii) small critical exponents 
r associated with the power-law fragment-size distribution, and (iii) the typical pattern of cracks. 
The proposed computer simulations do not require numerical solutions of the Newton's equations of 
motion, nor several additional assumptions normally used in discrete element models. The model is 
also able to predict some physical aspects which could be tested in new experiments of fragmentation 
of brittle systems. 

PACS numbers: 46.50.+a, 62.20.Mk, 91.60.Ba 



I. INTRODUCTION 

The scientific and technologically important subject of failure, fracture, and fragmentation in condensed matter 
physics has a rich phenomenology and depends on a broad range of physical and chemical factors such as material 
composition, impurities, defects, internal and phase boundaries, and temperature, among others. Because of the 
complex interactions involving these factors, in most cases it has been difficult to develop a consistent understanding 
of fragmentation processes on a fundamental level. To complicate matters further, statistical complexities arise in 

^ : n 

experiments on fragmentation 1]. For instance, the experiments are invariably statistically inhomogeneous; i. e. the 
fragment size varies as a function of position within the fragmented object because of the spatial variation of forces 
and energies causing fragmentation. The concept of scaling is important in a number of fragmentation processes 2|, 
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and critical exponents can be precisely defined in several cases, depending |3|, or not J] on the shape of the object. 
More recently, the fragmentation of low-dimensional brittle materials embedded in a higher dimensional space has 
been studied by several groups [5|, [fl LD, Is |8 , 
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In this work we use the random capability growth (RCG) model, a simple computer simulation whose rules are 
inspired by population genetics and which was recently introduced to study the evolution of the linguistic diversity 



on the Earth 



12, 



13], to describe the statistical details observed in experiments dealing with brittle platelike objects 
fragmented under several types of impact P, 0] . In impulsive fragmentation, the RCG model includes the effects of 
inhomogeneity and randomness of the brittle material, as well as the effect of the diffusion of the stress pulse causing 
the fragmentation. Furthermore, the RCG model predicts new physical aspects that could be tested in experiments. 



II. THE MODEL 

The RCG model is defined on a two-dimensional lattice composed of A — L 2 sites with periodic boundary conditions, 
which simulates a two-dimensional plate whose height h « L. We assume that due to the impact there is a stress 
pulse that diffuses or disperses in the system starting from the impact point [14]. During this diffusion process some 
connections between lattice sites are broken depending on their particular mechanical strength. As we are assuming a 
disordered non-crystalline medium, this mechanical strength is a randomly distributed variable on the lattice. Thus, 
to each lattice site i we ascribe a given capability Ci whose value we estimate from a uniform distribution in the 
interval 0-1, and so Ci refers to the facility to accomodate at site i an excess of stress from the original impact. 

In the first step of the dynamics, we randomly choose one site where the excess of stress due to the impact is firstly 
injected. This site will belong to the fragment labeled by number 1. In the second step, one of the four nearest 
neighbors of this site will be chosen to receive the stress pulse with probability proportional to Ci (see Figure la). 
We consider that there is a finite probability p to occur a mutation which will break a single connection between two 
sites and which will contribute to the formation of a new fragment at a future time. In the context of population 
genetics mutations are the mechanisms responsible for generating diversity. In the present case, if a mutation occurs 
the chosen neighbor site is labeled by number 2 and the nucleation of a new fragment begins (Figure lb). Otherwise, 
it is labeled by number 1 and the size of the initial fragment is increased. The mutation probability is given by 
p = a/S, where a is a constant in the interval 0-1, and S is the fitness of the fragment, defined as the sum of all the 
capabilities of the sites belonging to that specific fragment. Therefore, the initial fitness of the first (impact) site is 
the capability of the initial site. This rule for the mutation probability was inspired by population genetics where the 
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FIG. 1: In a lattice composed by A — 7 2 sites we show: (a) Impact site (labeled by number 1) and its four nearest neighbors 
(time t = 1); (b) The occurrence of a mutation begining the nucleation of a new(future) fragment labeled by the integer 2 (time 
t = 2); (c) Cluster whose sites were visited by the stress pulse and its boundary at time t = 11; (d) Increasing of the fragment 
labeled by the integer 1 (time t = 12); (e) The begining of the nucleation of a new fragment which is labeled by the integer 3 
(time t = 12). 

most adapted organisms are less likely to mutate than poorly adapted organisms [15]. The increase of the parameter 
a introduces an increment in the mutation probability, that is, a is a direct measure of the impact energy. 

In the subsequent steps, we check the sites which are on the boundary of the cluster whose sites were visited by 
the stress pulse, and we choose one of those sites according to their capabilities (Figure lc). Those with higher 
capabilities have a higher likelihood to be visited by the stress pulse. We then choose the fragment to which this site 
will belong among their neighboring fragments. The fragments with higher fitness have a higher chance to increase in 
size (Figure Id). As before, the mutation probability is inversely proportional to the fitness of the chose fragment and 
if a mutation occurs the site is labeled by number 3 (Figure le). This process will continue until all sites be visited 
by the pulse. At this point, we verify the total number of fragments N. Each selection of a site corresponds to an 



4 




FIG. 2: (a) Snapshot of a particular realization of the dynamics after all sites be visited by the stress pulse. The lattice size is 
L = 20 and a — 0.6. (b) Typical distribution of fragments observed in a glazed-tile after impact with a conic steel projectile 
falling in free fall. See text for detail. 
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FIG. 3: Number of fragments iV as a function of the area A for a = 0.6. The exponents are v = v+ = 0.69 ± 0.01 for 
4 < A < 7000, and u = v- — 0.22 ± 0.01 for 7000 < A < 250, 000. In the inset we present exponent vasa function of a for (a) 
small and intermediate areas and (b) large areas. 



increment of one time unit. 

With the passage of the time, the probability increases to find fragments with increasing fitness, thus when the 
distance from the (initial) site of impact grows, the size of the fragments increases as well, as physically expected. 
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FIG. 4: Number of fragments with area greater than a, n(> a), as a function of a, for a — 0.6, and L = 500. n(> a) ~ a _s 
with B = B + = 0.60 ± 0.01 for 1 < a < 20 and B = S_ = 0.27 ± 0.01 for 45 < a < 12, 000. In the inset we present the 
exponents B- (bottom) and B+ (top) as a function of a. 

Figure 2a shows a particular configuration of fragmentation on a lattice with 20 2 sites. For comparison, Figure 2b 
shows a typical distribution of fragments observed in a glazed-tile of size 150 x 150mm 2 and thickness 5mm, after 
impact with a conic steel projectile of 238g of mass falling in free fall from a height of 5m 



16] 



III. RESULTS AND DISCUSSION 



Due to the effect of increasing fitness for large times, discussed in the previous section, if we fix the impact energy 
(a fixed), and for plates of increasing areas, we observe that the number of fragments per area, N/A, is a decreasing 
function of A. This is exemplified in Figure 3 which shows a composite power-law for N(A) for a simulation with 
a = 0.60, and average on 200 similar experiments for L < 300, and on 20 similar experiments for L = 400, and 500: 



N~A", (1) 

with the exponent v assuming the value v+ = 0.69 ± 0.01 for small and intermediate areas, and v_ = 0.22 ± 0.01 for 
large areas. The dependence of the exponent v with the parameter a is shown in the inset of Figure 3. As a increases, 
v increases towards the maximum value v = 1 in an approximately linear way, and the composite scaling reduces to 
a single scaling when v = 1. Unfortunately fragmentation experiments in general do not examine the dependence of 
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FIG. 6: Number of fragments with area greater than a, n(> a) as a function of a for a = 0.6 and L = 500 for fmax = 1000 (□) 
and f m ax = 10 (*). When f max = 1000 we estimate B+ = 0.61 ±0.01 for 1< a < 30 and B- = 0.45 ±0.01 for 30 < a < 10, 000. 
When fmax = 10, we have B+ = 0.71 ± 0.01 for 1 < a < 30. 

the total number of fragments with the area of the plates. Thus after the experimental measurement of N(A), we 
could compare the data with our prediction in equation 1, and Figure 3. 

If we work with a scaling distribution of fragments of the type n(a) = n(l)a~ T (where n(l) is the number of fragments 
of minimun (unit) area), or an accumulated distribution of the type n(> a) ~ a~ B (B = r — 1), both defined in the 
interval 1 < a < a max on a 2D system of area A = L 2 , we have the sum rule A = J n(a) ■ a ■ da = 
The last result reduces to A w [n(l)/(2 — T)]a^J x , (for r < 2, and a ma x >> 1, as commonly observed in 2D 
real fragmentation processes for not too large impact energy). Furthermore, based on many experimental results 
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FIG. 7: Evolution of N(t) for a = 0.6 and L = 500. In this curve we have N(t) ~ t° 
JV(t) ~ f 08 ± 01 for 24, 000 < t < 250, 000. 



for 60 < t < 2, 100 and 



it is valid to assume that the largest fragment size a max is proportional to the total area A of the system (also 
for not too large impact energy); in this case we get from the last expression n(l) ~ A 1 " -1 . Thus, we find that 
N = J n(s)ds « [n(l)/(r — 1)] ~ A T_1 = A s , provided r > 1; i.e. if the scaling (1) is valid, the expected differential 
distribution of fragments n(a) tends to scale as a ~ v or, equivalently n(> a) ~ aT v , with v = B. 

The most frequently studied statistical function in fragmentation experiments is the average number of fragments 
of size larger than s, n(> s), or area larger than a, n(> a), for breakup of plates. The last statistical function for 
our fragmentation model is shown in Figure 4 for a — 0.60, and L — 500 (10 similar experiments). In our model, 
n(> a) exhibits in general a composite scaling law characterized by two different critical exponents B+, and _B_ valid 
respectively for small and intermediate area fragments. Thus, the expected differential distribution of fragments, n(a), 
which is the derivative of n(> a) with respect to a, scales as n(a) ~ a~ T , with r assuming the values t + = B + + 1, and 
t_ = B-+1, for small and intermediate areas, respectively. With the power law fit in Figure 4 we find B + = 0.60±0.01 
(r+ = 1.60 ±0.01), and B_ = 0.27 ±0.01 (r_ = 1.27 ±0.01). It can be noticed that the exponents B (or r) obtained 
from Figure 4 are in agreement with the relation B = v and with Figure 3. For the reader to develop a better insight 
into our model, we show B as a function of a in the inset of Figure 4. We notice from this inset that both B_ and 
B + present a nearly linear dependence with a, with a constant difference {B + — BJ). 

The behavior observed in Figure 4 for n(> a), with composite scaling laws, is in very good agreement with the 



experiments of Meibom and Balslev 



4[ of fragmentation of plates of dry clay falling onto a hard floor from a height 



of 2.0m (In 4| the experimental cumulative functions corresponding to n(> a) are shown, although the critical 



8 

exponents reported there, /3, are those of the differential distribution, i.e. the same as t, in the standard notation 
adopted in the present work.). In the experiments of Meibom and Balslev, the observed values for the exponents were 
t + = 1.57 ± 0.09, and r_ = 1.19 ± 0.08, irrespective the aspect ratio of the plates, h/L, in the interval 0.015 to 0.133. 
It is interesting to observe that the exponent t_ = 1.27 ± 0.01 found in Figure 4 is only 5% off from the exponent 
t = 1.35 ± 0.02 obtained in recent experiments of impact fragmentation of egg shells [5]. Furthermore, for a = 0.10, 
the accumulated distribution is shown in Figure 5. In this case, we observe the power law n(> a) /a ~ a -5-1 = &~ r , 
with t = 1.08 ± 0.01 exactly as observed by Oddershede et al [3J with the fragmentation of gypsum disks of 320mm 
in diameter and 5mm thick. In our Figure 5 the scaling persists over four decades in a, although in Figure 2 of 3| 
the scaling occurs only over two decades due to the difficulty to detect experimentally those fragments with mass in 
the range 5 x 10~ 5 g to 5 x 10 _3 g. 

In order to understand the role of the brittleness of the material on the critical exponents, we assumed that the 
fitness of each fragment is bounded by a given maximum value, f max , which is randomly chosen from an uniform 
distribution. In this way, the smaller f max , the more brittle is the material. In Figure 6, we present the average 
number of fragments of area larger than a, n[> a), for a = 0.6 and L — 500, when f max = 1000 and f m a X = 10. 
When we compare these results with those shown in Figure 4 (where there is no limit to f m a X ), we notice that the 
exponent B + is roughly the same when f max = 1000, and a slight deviation of 10% is obtained for f max —W- This 
seems to indicate that B + is insensitive to the brittleness of the material. On the other hand, we see that i?_ varies 
with the brittleness. For instance, when f max — 1000 we estimate = 0.45 ± 0.01, which represents a deviation of 
50% in comparison to the non-bounded case (or a deviation of 14% in t_ = 1 + BJ), and for very small values of 
fma X , the scaling region described by £?_ does not exist. 

Finally, Figure 7 shows the time evolution of the number of fragments, N(t), for the same values of a and L used 
in Figure 4 (average of 10 experiments). The time scaling appearing in Figure 7 could be tested in new experiments 
of fragmentation of brittle platelike objects using present-day high-speed high-resolution digital imaging techniques. 

IV. SUMMARY AND CONCLUSION 



It can be observed that the RCG model is useful in obtaining not only the composite scaling laws observed in the 
fragmentation of brittle materials [4(, but also the small exponents r (close to 1) 3J, |4|, |5(, and the typical pattern of 
cracks observed experimentally (Figure 2). These achievements of the RCG model are obtained without resorting to 
out-of-plane bending modes and discrete element models where the dynamics of the breakup process are determined 



by molecular dynamics simulations, as in references [5| and [8|. Thus it is pertinent to ask: Why indeed should the 
simple RCG model exhibit all these aspects observed in experiments of fragmentation with brittle materials? This is 
not an easy question. A possible answer, though, is suggested just by the two ingredients of the model: It could be 
said that (i) diffusion (in a disordered medium described by random capabilities) of a stress pulse due to the impact, 
and (ii) breakup probability stochastically controlled by the fitness function (as described in the fourth paragraph) 
are very important physical ingredients to be considered in the understanding of the basic statistical aspects of low- 
dimensional brittle fracture. It is interesting to note that of these two ingredients inspired by population genetics 
simulations, the ingredient (i) is physically plausible in a number of fragmentation processes, while the ingredient 
(ii) is more subtle, and deserves a further detailed investigation which is outside the scope of the present paper. 
Moreover, it is opportune to stress that the understanding of the physical origin of nonequilibrium structures like 
the crack patterns obtained in impulsive fragmentation remains one of the most challenging problems in statistical 
mechanics. 

In conclusion, we have used a simple one-parameter fragmentation model which describes quite well the rich 
phenomenology observed in the fragmentation of brittle platelike systems, and in particular the existence of composite 
scaling laws. The model predicts new scaling laws, and we urge that more experiments be performed to investigate 
the laws proposed here. 

We thank I. R. Tsang by the Figure 2a, V. P. Brito for provide us with his experimental data showed in Figure 
2b and M. L. Sundheimer for proofreading the manuscript. V. M. de Oliveira and M. A. F. Gomes are supported by 
Conselho Nacional de Desenvolvimento Cientffico e Tecnologico and Programa de Niicleos de Excelencia (Brazilian 
Agencies) . 
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